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The interaction of the flow through the labyrinth seals with the shaft of the rotor can have 
an effect on the stability of turbomachines. Thus, the excited forces, so-called cross forces or 
non— conservative forces, arise, which act perpendicular to the rotor eccentricity. This effect 
is caused by an unsymmetrical pressure distribution within the labyrinth cavities. 

Experimental investigations are carried out for different types of labyrinth geometries: 
two staggered labyrinths with teeth on the stator and grooved rotor as well as a full and a 
convergent stepped labyrinth. These labyrinths can be found on the tip shrouding of bladings 
in steam or gas turbines. The following parameters are varied in the test-facility: geometry 
of the labyrinth seals (number of cavities, inlet region), shaft rotation, pressure difference on 
the seal, entry-swirl and eccentricity of the rotor. 

The results are presented for stiffness coefficients of the labyrinth seals, leakage flow and 
circumferential flow in each cavity which was measured with special probes. Generally, the 
inlet swirl has the greatest influence on the coefficients of the seals. The experimental results 
are compared with theoretical results and are in good agreement. 


*Work under this study of labyrinth seals has been sponsored partially by the 
Forschungsvereinigung Verbrennungskraftmaschinen e.V. (Research Association 
for Internal Combustion Engines), Frankfurt/Main, Federal Republic of Germany. 
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1 Nomenclature 


A 

B 
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u = c/Rru 

U 

z 

/?* 

e - e/CR 


Area 

Height of labyrinth seal strip 
Stiffness-term 

Circumferential velocity of flow 

Inlet swirl 

Damping-term 

Diameter 

Flow characteristic 
Dimensioned eccentricity 
Force 

Reference force 

Cross-sectional area of control volume 

Cross-sectional area for centered shaft-position 

Value of parameter of Kq 

Local radial clearance 

Hydraulic diameter 

Complex unity 

Factor of area change 

Cross force spring-coefficient 

Restoring force spring-coefficient 

Inertia-term 

Number of cavities 

Coefficients for friction factor 

Leakage flow 

Leakage flow per circumferential length 
Leakage flow for centered shaft-position 
Pressure 

Pressure before seal 
Position of probe 
Cross force 
Restoring force 
Gas constant 
Reynolds number 
Shaft radius 
Nominal clearance 
Temperature 
Time 

Pitch of seal strips 

Rotor peripheral velocity 

Dimensionless circumferential velocity of flow 

Length on which shear stress acts 

Number of seal strips 

Flow coefficient 

Eccentricity ratio, perturbation parameter 
Flow coefficient 
Kinematic viscosity 


no 



p Density of fluid 
7 r Pressure ratio (p a /po) 

ip Peripheral angle 
r Dimensionless time 
a Shear stress 
lo Shaft angular velocity 
Q, Shaft precessional velocity 
subscripts 

0 zeroth-order component 

1 first-order component 
c cosine 

1 i-th chamber value 
Q cross 

R rotor, restoring 
r radial-direction 
S stator 
s sine 

u circumferential direction 
w shaft 
y y-direction 

2 z-direction 
- mean-value 

* dimensionless value 
Not presented notations are declared in text. 

2 Introduction 

In order to carry out a thorough stability analysis of turbomachine rotors, it is essential to know all the 
forces causing oscillations. These forces, acting on the rotating subassembly units of the rotor shaft, can 
effect the stability behaviour in such a way that the operation is seriously disturbed or even impossible. 
The major causes for rotor instability are mostly exciting forces resulting from fluid-mechanisms, which 
are found in journal bearings and in labyrinth seals. The interaction between the flowing medium and the 
rotor shaft cause a transmission of energy to the flexural oscillation system. That leads to rotor-instability, 
if damping is not sufficient. 

In order to be able to judge the stability qualities in advance, the engineer has to carry out oscillation 
analyses. For that he needs information about the dynamic characteristics of the seals which are charac- 
terized by stiffness, damping and inertia terms. When these properties are known the fluid-mechanisms 
can be described in a linear equation of motion, i.e., equation (1). 



These quantities must be known in order to design the dynamics of turbomachines properly. 

Experimental results of labyrinth seal forces have been presented by Benckert and Wachter [1], [2]. 
They investigated several multistage types of different labyrinth seals and found out that for long seals 
the shaft rotation has a great influence. 

Childs and Scharrer [3] surveyed experimentally the dynamic characteristics of long straight through 
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labyrinths. On this test facility, the stiffness behaviour as well as the damping characteristics could be 
determined, since the rotor carries out an oscillation motion in addition to its rotational motion. 

Leong and Brown [4] surveyed long straight through labyrinths with different cavity geometries. The 
seal strips were attached to the stator. Their results show, that with increasing peripheral velocity of the 
rotor, the cross-coupled stiffness coefficients decrease and can become even negative. 

Rajakumar and Sisto [5] and Serkov [6] investigated the flow forces for short straight through labyrinths. 
They found out, that the inlet swirl compared to the shaft rotation has a great influence on the cross forces. 

In reviewing previous experimental investigations, there is a clear need for testing of short prevalent 
labyrinth seals, which can be found on the tip shrouding of bladings in steam or gas turbines. 


3 Test— Facility 

The investigations are carried out on a labyrinth test-facility at the ’Institut fur Thermische Stromungs- 
maschinen und Maschinenlaboratorium’ of Stuttgart University. The test-facility was originally set up in 
[1], [2] and has now been modified to fit the changed requirements. Using this setup the inlet flow is the 
same as in the shroud ring seals of a real turbine stage (e.g., overlapping of guide- and rotor-bladings, 
radial streamline variation). 

The cutaway (figure 16) shows one of the shroud ring seals with grooved rotor and teeth on the stator 
to be examined. The flow force in the labyrinth seals is always determined by integration of the measured 
pressure in axial direction andlalong the perimeter.) Inlthelaxial|directionieach|cavity|isllookedlat individually 
and then all cavities are summed up. The reduction of these results to components gives the cross and 
the restoring force Q and R. The influence of the following parameters have been surveyed: 

• Inlet state before the first seal tip tc,c u q. It depends on the pressure ratio being at the seal and the 
inlet swirl of the flow. A nozzle ring in conjunction with a vacuum assembly can produce different 
swirl components c u o- The measuring of the resulting circumferential component of the labyrinth 
flow is made with three-hole- cylinder probes specially developed for this purpose (see figure 1). 
Since the probes are also used in the area of the labyrinth strips and near the rotor surface, the wall 
effect on the calibration factors, defined in figure 1, must be tested and taken into consideration. 
The test medium is air which expands to the ambient conditions. 

• Eccentricity of the rotor e. The eccentricity is varied from negative over the centered position to 
positive values. 

• Peripheral velocity of rotor u w . The rotational rotor speed is increased from u w = Om/s upto over 
u w = 100 m/s. The influence of the direction of rotor rotation is likewise being tested. 

• Geometry of the seal. Measurements are made with 15 different labyrinth configurations. They can 
be subdivided into three different types of labyrinths; grooved, interlocking and stepped labyrinths 
(figures 18). Besides the number of cavities, the inflow contour of the labyrinth was also varied. 

Figure 2 shows a summary of variations of the test parameters. 
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Figure 1: Three-hole-cylinder probe 


Variation of Test Parameters 
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Figure 2: Variation of test parameters 
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4 Analysis of test data 

The results are presented as follows: 

» Variation of the flow forces in the individual cavity as dimensionless cross and restoring force in each 
cavity. The restoring force R* acts against the rotor eccentricity e (centering effect) and the cross 
force Q* acts perpendicular to e. 

• Variation of the forces in the entire labyrinth as dimensionless cross force spring coefficient Kq 
dependent on the test parameters and the flow characteristic Eq. The quantity Eq is the ratio of 
the inlet swirl energy and axial pressure drop. 

* Leakage characteristic 

The definitions for the calculation of the forces and the stiffness coefficients are summarized in figure 17. 


5 Experimental Results 

Forthe|different|labyrinth|typleslof M- and A-configuration, the figures 19, 20 show the pressure distribution 
and the resulting rotor loads, i.e. cross and restoring force Q* and R*. The pressure distribution p*(<p) 
is related to the static pressure difference. For three eccentricities (e = 0 and e ± 0,66), the pressure 
distributions depending on the individual labyrinth cavities and its circumferential angle (p are depicted. 
The flow conditions r,c u0 ,u w are noted 1 - in the figures. 

For centered position (e = 0), the pressure distributions along the perimeter of the rotor are almost 
constant in the individual cavities. Thus, the resulting forces Q* and R* are nearly zero. But looking at 
thepressure|drop in the axial direction|shows,|that|the| pressure drop depends on the geometry of the cavity 
(transitions of the staggered grooves). Figure 3 shows the resulting flow coefficient //,■ as a function of the 
cavity number i for the Ml to M4-configuration. 



”» / PqVq 
A <\ Pi - pf +1 



Figure 3: Flow coefficient m of Ml-M4-configurations 
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For positive and negative eccentricity (e± 0, 66), the variations of pressure | in 1 the circumferential 
direction are almost identical. This correspondence for ±e, (which must comeaboutifor a perfect labyrinth 
geometry and an exactly adjusted rotor, is even more evident for Q* and R* . 

The results show, that it is not the first cavity which provides the highest contribution of cross force. 
In the front part of the labyrinth, the cross force rises to its maximum and drops then continuously. 
For the M2-configuration the maximum is located in the 2. cavity (see figures 19, 21). This pattern is 
independent of inlet swirl, eccentricity and pressure rate. It was detected for all M-configurations. But 
the increase of the cross force in the front part depends on the inflow contour of the labyrinth. Relatively 
high gradients dQ*/di were measured for the inflow contour with rotor groove and opposite long seal tip 
(M2 to M4-configuration). In figure 4 this is depicted for the M3 and M4 -configuration. 




X 

e=+0.66 

Cu O =110m/fe 



V 

£= 0 

u w = 0 m/s 



+ 

£=-0.66 

7t=0.57 


. 0.2- 
o 0.18 + 
0 0.16- 



X 

M4, 

O 0. 14-; 

p 0.12- 


X 

+ X 


h* 0.1- 


+ 



*8 0.08- 
M 0.06- 

t/5 0.04- 
O 0.02- 


¥ 





V v 

ro 

II 

E 

O 


1 1 

i i 



ges 1 

2 3 




Cavity 




Figure 4: Cross force Q* of M3- and M4-configuration 


The increase of the cross force dQ*/di is smaller for a inflow contour with rotor stage and opposite 
short seal tip (Ml-configuration, see figure 5). 
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Figure 5: Cross force Q* of Ml-configuration 
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Apparently, the changed inflow geometry brings about a different inflow condition at the labyrinth. 
An extension of the space in front of the first seal tip results in a larger increase of the cross force in the 
front part of the labyrinth. 

Like the experiments, analytical investigations show different cross force characteristics for each cavity 
when varying the boundary conditions. The boundary conditions are: the pressure distribution and the 
circumferential velocity distribution of the labyrinth flow along the perimeter in front of the first seal tip, 
and the flow coefficient jig. Furthermore, the boundary conditions at the exit of the labyrinth may also 
have, an influence. 

The influence of the inflow geometry on the cross force characteristics could be proved for the A-full- 
labyrinth as well (see fig. 6). 




Figure 6: Cross force Q* of A2- and A 3-configuration 

For the M2-configuration, fig. 21 shows the influence of the inlet swirl c u o and the eccentricity e on 
the cross force- and restoring force characteristics. 

In contrast to Q *, the restoring force is alternating along the cavities. Hence the resultant restoring 
force R *, recorded at the cavity number i = ges, is almost zero (see fig. 7). 



Figure 7: Restoring force R* of Al-configuration 
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This characteristic is even more evident for the A-full-labyrinth. Pressure maxima and -minima 
respectively dependent on the cavity number occur in the area of the narrowest gap (p = 180°). Similar to 
the Lomakin-effect [7] , the losses due to friction and inflow respectively at each seal strip varying along the 
circumferential direction are responsible for that. That means, besides the dependence on the arrangement 
of the seal strips, the flow coefficient [m is also a function of the peripheral angle <p. 

The alternating run of the restoring force R* is independent from the peripheral velocity of the rotor. 
This is shown in figure 7 for the A 1-configuration with 7 cavities and c u o = 0 m/s. 

For the Ml-(m=8) and M4-configuration (m=3) the figure 8 illustrates Q* versus cavity number i for 
different u w and superposed c u o- 
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Figure 8: Cross force Q* of Ml- and M4-configuration 

If the shaft rotates in the opposite direction of the inlet swirl, the value of u w becomes negative. By 
definition the inlet swirl is positive in the direction of the peripheral angle (p. Changing the direction of 
u w different situations in turbines and compressors can be simulated. Figure 8 shows that the cross force 
for u w = —95 m/s is higher compared to u w = 0 mjs. 

This effect can be explained by the ratio of decelerating to accelerating friction areas at rotor and 
stator. It has an effect on the distribution of the circumferential component of the labyrinth flow along 
the seal. The variation of c u along the cavities shows that the amount of cross force doesn’t depend on 
the amount of circumferential velocity of the labyrinth flow, but on the change of the circumferential flow 
between the cavities. For negative u w the swirl changes are maximal (see fig. 9) and consequently the 
cross force as well. 

With the help of flow-measuring probes it was ascertained that the circumferential component of the 
labyrinth flow is nearly not radial-changing within the cavities. Fig. 10 illustrates c u versus the cavity 
number i for centered shaft position. Only in the vicinity of the stator surface the c u -components are not 
constant. This is because of the effect of the boundary layer. Hence the circumferential velocity of the 
labyrinth flow in the cavities can be approximated by a mean circumferential velocity. 


117 





M2, n=0.47, e=0 


x iv»+10Q m/« * iV- 0 m/« 


♦ u^-100 m/« PI - 0.7 

no- 

+ v X 

90 — 

x 


* X 

70- 

* V X 

r 5T 

\ 

JL 50- 

+ * X 

3 


30- 

+ 

10- 


—10 — 

1 + 




Cavity 



Figure 9: Circumferential velocity c u with superposed rotor speed 
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Figure 10: Circumferential velocity c u of M 4-configuration 

Even for eccentric rotor position the c u ~components are almost constant in radial direction. Further- 
more, fig. 25 shows for the M 4-configuration the| variationiof.c u j in, the|peripherall direct ion| in each cavity. An 
approximation of c u in terms 

c Uj i — ^um t i "f £ cos (p -|- sin<p) (2) 

confirms the assumptions of equations (9) and (11), see figure 25. 

The relationship between the total cross force Q* and the eccentricity e is linear for all configurations 
as the plots in figure 11 show. 

The characteristic field of the spring coefficient Kq as a function of ir, c u o and u w is shown in figure 
22 and 23. These characteristic fields enable a separate investigation of the varied test parameters and its 
effects on the spring coefficient. It shows that the inlet swirl has the greatest influence on Kq. 
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Figure 11: Cross force Q* versus eccentricity e of Ml -configuration 

Introducing the flow-characteristic Eq, the results can be transferred to seals in real machines. Fig. 
1*2 shows for the M-configurations the obtained spring coefficients and the approximated curves. The 
approximating function is given by the nonlinear equation Kq = gi+giE^ 93 . The coefficient gi corresponds 
to spring coefficient Kq for axial inflow ( Eq = 0). 




Figure 12: Kq versus Eq of M-configuration 
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The comparison of the four M-configurations shows that the parameters < 71 , 2,3 depend on the number 
of cavities (see table 1). 

In addition, table 1 shows the flow coefficient f3 * for Eq = 0. The relationship between the leakage 
flow m and the flow coefficient /3* is: 


771 



( Pa-Pi 

Po/po 


Confi- 

gurat. 

Uw 

[m/s] 

5i 

52 

53 

(3* 

Ml 

0 

-0,0031 

0,3551 

0,7336 

1,9 

Ml 

47 

-0,0092 

0,3809 

0,7687 

1,86 

Ml 

94 

-0,0115 

0,4174 

0,8252 

1,79 

M2 

0 

-0,0023 

0,3533 

0,6591 

1,84 

M2 

47 

-0,0081 

0,3210 

0,6489 

1,82 

M2 

94 

-0,018 

0,298 

0,60 

1,76 

M3 

0 

-0,0016 

0,4726 

0,6485 

1,78 

M3 

47 

-0,0139 

0,4843 

0,6710 

1,76 

M3 

94 

-0,0327 

0,5091 

0,6636 

1,72 

M4 

0 

0,0008 

0,7586 

0,8059 

1,68 

M4 

47 

-0,0109 

0,5063 

0,6472 

1,67 

M4 

94 

-0,0326 

0,4826 

0,6004 

1,64 

Al 

0 

0,0030 

0,7316 

0,8132 

1,56 

A1 

47 

0,0004 

0,6954 

0,7743 

1,55 

Al 

94 

0,0062 

0,7616 

0,8274 

1,53 

A2 

0 

0,0020 

0,7747 

0,8978 

1,55 

A2 

47 

0,0010 

0,7132 

0,8362 

1,54 

A2 

94 

-0,0007 

0,6859 

0,8046 

1,52 

A3 

0 

0,0044 

1,0057 

0,8812 

1,45 

A3 

47 

-0,0013 

0,7879 

0,7544 

1,44 

A3 

94 

-0,0058 

0,7110 

0,6955 

1,42 


Table 1: Values of the parameters < 71 , 2,3 Kq and /?* of M- and A-Iabyrinth 

Furthermore, the characteristic fields of the restoring force spring coefficient as a function of the 
test parameters are depicted in the fig. 24. 
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6 Analysis 

In order to compare the experimental with the theoretical results an analysis was developed. It bases on 
Iwatsubo’s [8] and Child’s [9] basic equations for compressible flow in a labyrinth seal. The continuity and 
momentum equation will be derived for a single cavity control volume as shown in figure 13. 
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Figure 13: Cavity control volume 


The governing equations are 

• Continuity equation: 

0 \ 3 

Ft + (HA*) + ~ = 0 

• Momentum equation: 
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dt 
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• Leakage equation: 

— piHi 

with leakage mass flow per circumferential length 

• Friction equation: 

1 2 . 1 2 f\Ci\h di \ msi . 

&t,Si = 2 P' c i*Si = 2 Pi c i n Si ( — ~ — ) sign (c t -) 
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with hd as hydraulic diameter. 



( 3 ) 


( 4 ) 
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Ideal gas law: 


— = RT 
Pi 


® Normalization of the variables: 
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Rru Po Po 'Mq 

The solution procedure is similar to the perturbation analysis used by Childs [10]: 

• Linearized zeroth and first-order perturbation equations describe the pressure and velocity distribu- 
tions for the centered (index ’O’) and eccentered (index T’) shaft position. The following| pertun- 
batioh variables are defined: 

U'l — ^0 ,» T Pi — P0,i T fPl,j fi — fo,i T tfl,i 

Hi = Cr + eHi t i pi = po t i + epi'i (9) 

where e = e/C r is the eccentricity ratio. 

• The zeroth-order solution for the circumferential velocity (zeroth-order momentum equation) is 

RTM 0 
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( 10 ) 


The partial differential equations (first-order solution) will be reduced to linear complex equations by 
introducing the following assumptions for the perturbation variables: 


• Elimination of circumferential angle < p: 

Pl,i = 
hl,i — 

• Complex variables: 


Pi,i,c cos <p + sin ip Ui t i = Ui t i iC cos ip + u it i iS sin ip 

H u 


Cr 


= -yi cos ip - Zi sin p> /*!,*' = p ltiiC cos ip + Mi,»> sin <p 


( 11 ) 


P.1 i — Pl,i,c "b jPl,i,s — ^l,t,c T j u l,i,s 
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( 12 ) 


• Elimination of time t with harmonic seal motion (circular orbit with frequency ft around the cente- 
red position, see fig. 14) : 


u . _ l .Jttt 
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and corresponding response for the other perturbations variables: 
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( 13 ) 


(14) 


yields to the complex continuity and momentum equation: 
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Complex continuity equation: 
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(16) 


• With equations (15), (16) a system equation for all cavities is formed and solved by linear equation 
algorithms. As a result the values of the complex amplitude p^ . and u ljt - (p i = p 1)!C 4- jpi,i s ) are 
found. 

• The flow force for the seal (see equation (1)) is found by integrating the first-order pressure pertur- 
bation along and around the shaft of the rotor. For the i-th cavity it is 


F y< i(t ) = -P 0 Rr£ 7r J (p\ t i c cosClt — pi t i s sin ftt) dx 
0 

F z ,i(t) = -PqRr£kL (pi,i, c sin ft! + p lti>s cos fti) 

Introducing a rotating coordinate-system (see figure 14) the rotating forces are 

N 
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Figure 14: Harmonic seal motion and transformation of the coordinate system 




With the assumed orbit about the centered position the structure of the matrices of equation (1) 
becomes 


Cyy — C-zz i C : 


J/2 


— C. 


yz) 


* yy 


d'ZZ ) dy Z dyz , 


(19) 


• The searched coefficients can be found by a least square fitting for different frequencies of ft: 


F r — cCf i^yy “1“ dyz'H Tflyy^l ^ 

Fy, = -eC r (c yz - dyyQ - my Z Q. 2 ^) (20) 


7 Comparison of the Results 

Comparisons of theoretical and experimental results of the exciting cross-coupled stiffness coefficient c yz 
respectively spring coefficients Kq are made for the staggered labyrinth, configuration M4. The correlation 
is as follows: 

_ c yzCR 

Q Ap st R R mL 

The results are presented for Kq as a function of the flow characteristic Eq. The friction parameters (see 

equation (6)) for the evaluations are set to 

na = 0,079 m,Ri = — 0,25 nsi = 0,15 msi = — 0,25. 

The other input data like flow coefficient m and boundary conditions are derived from experimental 
results, see figure 3. 

The results are in good agreement with the spring coefficients Kq (see figure 15). 

Furthermore, an identification algorithm is currently under development. Its aim is to determine the 
friction factors n R -si and rriR-si for each cavity based on the zeroth-order solution for the circumferential 
velocity (see quation (10)) of the staggered labyrinth . The first results (see figure 23) confirm the 
assumptions of this bulk-flow-theory, equations (9), (11). 
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Figure 15: Comparison of theory and experiment 


8 Conclusions 

Experimental results have been presented for stiffness coefficients and leakage flow of different types of 
labyrinth geometries, such as grooved, interlocking and convergent stepped labyrinths. The influence of 
rotor speed, pressure ratio, inlet swirl and labyrinth geometry on the coefficients has been investigated. 

The results of this study support the following conclusions: 

• The cross force and restoring force are linear to eccentricity. Thus, the force-motion equation (1) is 
confirmed experimentally. 

• The inlet swirl has the greatest influence on the cross-coupled stiffness coefficient Kq . 

• The relationship between Kq and the flow characteristic Eq depends on the test parameters. 

• The cross force in each cavity depends on the change of circumferential velocity from cavity to cavity. 

• The restoring force shows an alternating behaviour along the seal. Thus, the total restoring force is 
nearly zero. 

• The results of the presented analysis are in good agreement with the experimental results. The 
experiments verify the assumptions of the analysis. 
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Figure 16: Testing facility of labyrinths 
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Figure 17: Definitions of test data analysis 
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Figure 20: Characteristic of p*(<p), Q*, R* of A1 
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Cross force Q- 



Figure 21: Q* versus c u0 , e of M2 



Figure 22: Characteristic field Kq of Ml 
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Figure 24: Characteristic field K% of Ml and M4 
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Figure 25: Circumferential velocity c u versus cavity i and peripheral angle ip of M4 
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